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H . Abstract 

<& . 

On a finite inhomogeneous graph model for complex network, we dehne the Ising model, 

which is a paradigm model in statistical mechanics. For the ferromagnetic Ising model, we 

calculate the Thermodynamic limit of pressure per particle. From our results, we compute 

other physical quantities such as the magnetization and susceptibility, and investigate the 

critical behaviour of this model. Our calculations use large deviation principles (developed 

Cln| recently) for suitably defined empirical neighbourhood measures on inhomogeneous random 

<-| graph. 
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,—l . 1. Introduction 

> 

— j Statistical ferromagnetic models or ferromagnetic ising models on Power-law random graphs such 

as the configuration model (CM), have received an interestingly increasing attention in the last few 

years , see e.g. Montanari et al. [2010] . Dommers et al. [2010] and the reference there in. For 

instance in Dembo and Montanari [2010] an explicit expression for the free-energy density and other 

important physical quantities such the internal energy in thermodynamic limits of n — > oo, where 

n is number of sites on the graph, were obtained if the degree distribution of the model has finite mean. 



* 



Their analysis makes use of the local neighbourhood properties of a randomly chosen edge which looks 
like a homogeneous tree where the offspring have size-biased degree distribution. The Ising model on 
a tree is simplier to analyse, and many intense scientific research work on it has been taking place. 
This is because the effective field of a vertex can be expressed in terms of that of its neighbours via 
the distribution recursion. The relation between the CM and a tree allows for the calculation of the 
internal energy and upon integrating over all possible values of the inverse temperature parameter we 
obtain the pressure. 

However, their technique cannot be applied to study the statistical mechanic on static inhomo- 
geneous random graphs (SIRG) or coloured random graphs, see Penman |1998J or Cannings and 
Penman [2003] which have the Erdos-Renyi graphs as a special case, because it is not locally-tree 
like with an asymptotic distribution. Some research, see e.g Agliari et al. [2010] . have use the Monte 
Carlo simulations to investigate the critical properties of the Ising model on Erdos-Renyi graphs. 
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In this work we use the LDP techniques developed in Biggins and Penman [2003], Doku-Amponsah 
and Moeters [2010J or Doku-Amponsah [2006J and furthered in Doku-Amponsah [201 lj to carry 
out Thermodynamic analysis of the Ising model defined on SIRG. To be specific, we obtain a 
thermodynamic limiting result of free-energy density via the random partition function. Using this 
results, we compute other physical quantities such as the, internal energy, the magnetization and the 
susceptibility of the system. 

Our main motivation for looking at these models is the study of the two-populations Ising model. 
The two-populations Ising model is relevant in many phenomena, ranging from the study of anistropic 
magnetic materials to social economic models, to biological systems. In the case of the latter , an 
interesting application concerns cooperativity effects in the context of bacterial chemotaxis, where 
lasting puzzle, namely how a small change in extrenal concentration of attractant or repellent can 
cause significant amplification in receptor signal. In theory, models onthe two-population Ising model, 
although simplified, are able to show that a coupled system of receptors has the capacity to greatly 
amplify signals. These values in limit, approache the critical coupling energy of an analogous Ising 
models. See, e.g. Agliari et al. [2010] and the reference there in 

2. Materials and Methods 

2.1 The Ising model. Let V be a fixed set of n sites, say V = {1, . . . , n} and E C 6 := {(-u, v) £ 
V x V : u < u}, where the formal ordering of links or bonds is introduced as a means to simply 
describe unordered bonds. Let X = {—1, +1} be the spin set and and denote by Q n {X) the set of all 
spinned graphs with spin set X and n sites. 

Given a symmetric function p n : R x R — > [0,1], a continuous function s : X — > R and a probability 
measure I on s(X) we may define the Inhomogeneous spinned random graph or simply spinned random 
graph X with n sites as follows: 

Assign to each site v £ V magnetic ising spin s(rf(v)) independently according to the spin law £. Given 
the spins, we connect any two sites u, v £ V, independently of everything else, with a bond probability 
p n [s(ji(u)),s(r)(v))] otherwise keep them disconnected. 

We always consider X = {(s(j](v)) : v £ V),E} under the joint law of graph and spin. We shall 
interpret X as spinned random graph and consider s(rj(v)) as the spin of the site v. 

On a spinned graph X, we define the ferromagnetic Ising Model by the following Boltzmann distribu- 
tions over X , 

HX{V) ~ Z x (/3,B) eX P\ 2^ fl(u) B(v) + Z^ VP J 

(u,v)eE uev 

where s(r)(u)) = \/]3B(r](u))r)(u), /3 > is the inverse temperature, B = {b(x) : x £ { — 1, l} n } is the 
vector external of magnetic fields, w(u) £ { — 1,1}, and Zx the random partition function is given by 

)/e{-i,i} v " (u,v)ge u&v 

The free-energy density or the pressure per particle is defined by 

$„(&£) = i log Z n (/3,B). 
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Our main concern in this paper is the study of the thermodynamic limiting behaviour of the free- 
energy, ie The behaviour of 

W,B):= hm 9 n {P,B) 

n— >oo 

at different level of temperature. Throughout the rest of the paper we assume that, for large n, X is 
sparse or near critical case i.e. The bond probabilities satisfy np n [s(x),s(y)] — > C[s(x),s(y)], for all 
x,y G { — 1, 1} and C: R x R — > [0, oo) is a continuous symmetric function, which is not identically 
equal to zero, and also satisfy the following energetic preference condition: 

(e/» _ i) ( C [-y/pB(-l), -v^S(-l)] - C[y/PB(1), y/pBQ.)]) = 2(B(-1) + B(l)). (2.1) 



2.2 Large deviation principles for spinned random graphs. In this subsection, we review some 
large deviation results of Doku et al. [2010] or Doku [2006] . and extend the joint LDP for empirical 
measures of coloured random graph to spinned random graphs. i.e. we assume a more general spin law 
£ : R — > [0, 1] with all its exponential moments finite and prove an LDP for this model in a topology 
generated by the total variation norm. 

To begin, we recall some useful definitions and notations from Doku et al. [2U10J or Doku [2006J. A 
rate function is a non-constant, lower semicontinuous function / from a polish space Ai into [0, oo], it 
is called good if the level sets {/(m) < x} are compact for every x G [0, oo). A functional M from the 
set of finite spinned graphs to Ai is said to satisfy a large deviation principle with rate function I if, 
for all Borel sets B C M, 

- inf Urn) < liminf ± logP„{M(X) G B\ 

meintB n->oo 

< limsup i logP n {M(X) G B) < - inf I(m) , 

n->oo meclB 

where X under P n is a spinned random graph with n vertices and int B and cl B refer to the interior, 
resp. closure, of the set B. 

And, for any finite or countable set y we denote by Ai(y) the space of probability measures, and 
by Ai(y) the space of finite measures on y, both endowed with the topology generated by the total 
variation norm. For uj G A4(y) we denote by ||w|| its total mass. Further, if £ G Ai(y) and u> <C £ we 
denote by 

#Mm =/<*] kg ($$) 

JM. 

the relative entropy of ui with respect to t. We set H(ui \\ £) = oo if uj ^ £. Finally, we denote by 
■M*(y x y) the subspace of symmetric measures in Ai(y x y). 

On each spinned graph X = ((s(rj(v)) : v G V),E) with n vertices, we define a probability measure, 
the empirical spin measure L 1 G At(Af), by 



L 1 i x \ ■= ~ ^2 § s(r,(v)) [x], for x G 



n 



and a symmetric finite measure, the empirical co-operative measure 1? G .M*(R x R), by 
L 2 [x,y] := - ^2 ( s {sMu), s( v (v)) + $(s( v (v), s(n(u)) )[x,y], for a,b <E X. 



n 

(u,v)£E 



The total mass ||L 2 || of the empirical pair measure is 2\E\/n. 
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Theorem 2.1. Suppose that X is a spinned random graph with spin law £ such that n" 1 \og£(n) — > — oo 
and link probabilities p„:Ixl4[0,l] satisfying np n [dx, dy] — > C[dx, dy], for all x, y £ R and some 
bounded symmetric measure C: R x R — > [0, oo). Then, as n — >■ oo, i/ie pair (L l ,L 2 ) satisfies a large 
deviation principle in M(R) x .M*(R x R) with good rate function 

1 



I(u, w) = /f (a; || /u) + 



2L 



i^ftu || Coo (S) w) + IICw <8> w| 



tu 



(2.2) 



where the measure Ceo ®w£ jM(RxR) is defined by Cuj <g> oj[dx, dy] = C(x, y)uj[dx]uj[dy]. 



2.3 Exponential Change- of- Measure. Given a bounded function /: X 
bounded function g : X x X — > R, we define the constant Ut by 



and a symmetric 



Uj = log 



„/[*! 



] £[dx], 



and the function /i„ : R x R — > R by 

fr n [x, y] = log (l - Pn [dx, dy] + p n [dx, dy]e 9[x ' y] ) 



(2.3) 



for a,b £ X. We use / and g to define (for sufficiently large n) a new coloured random graph as follows: 

• To the n labelled sites in V we assign spins from s(X) independently and identically according 
to the spin law £ defined by 

i[dx] = e f ' [x] - u f£[dx]. 

• Given any two sites u, v £ V, with u carrying spin a and v carrying spin b, we connect site u 
to site v with probability 



p n [dx,dy] 



1 - p n [dx, dy] +p n [dx, dy])e^ x ^ 



We denote the transformed law by P. We observe that £ is a probability measure and that P is 
absolutely continuous with respect to P as, for any spinned graph X = ((s(rj(v)) : v £ V),E), 



^ P ( v\ - TT felTM] IT Pn[ds(r,(u)),ds(ri(v))] 

dF y ] H Z[Mv(u)] 11 Pn[ds( V (u)),ds(v(v))} 11 l-pn[ds(r,(u)),ds(v(v))} 



TT l-p n [ds(r,(u)),ds(ri(v))] 

ueV (u,v)eE (u,v)£E 

n hds{y{u)} TT , ,., ..,. . . , ,, 
t[ds(r,(u)] 11 " < 

s_»uK„)-i- t rr e fl[*('/(«)i-M^'')). 

ueV (u,v)eE (u,v)€£ 

^ V {n{L\f-U I )+n{\L\~g)+n{\L 1 ®L\h n )-{\L\,h n )), 



p„[ds(ri(u)),ds(ri(v))] n-np n [ds(rj(u)),ds{r]{v))\ 
p„[ds(ri(u)),ds(ri(v))] n-np„{ds(ri(u)),ds(ri(v))} 

(u,v)eE (u,v)e£ 



n n-npn[ds(ri(u)),ds(ri(v))] 
n—np n [ds(r)(u)),ds(r)(v))] 



TJ e fHv(u)]-U f TT e g[s(v(u)),s(v(v))] TT e ^h n [s( V (n))Av(v))] 



(2.4) 



where 



L A = l^Z^X{u),X{n))- 



uev 



We write (/, w) := J M /[x]u;[dx] and 



(g,w) := 



#[x,y]tu[ofe,dy]. 
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Lemma 2.2 (Euler's lemma). If np n [dx , dy] —> C[dx,dy] for every x,y G R, then 

lim (l + apjd&.dyf)" = e aC[dx ' dy] , for all a,x,y el and /or a G R. (2.5) 

Proof. Observe that, for any e > and for large n we have 

(l + ggfefe )" < (l + ap n [<fo,dy])" < (l + aC ^ 1+£ ) n , 
by the pointwise convergence. Hence by the sandwich theorem and Euler's formula we get (|2.5p . ■ 

Lemma 2.3 (Exponential tightness). For every 9 > i/iere exists JVsN suc/i i/toi 

lim sup i log F\\E\>nN\ < -9. 

Proof. Let c > sup^ eR C[dx,dy] > 0. By a simple coupling argument we can define, for all 

sufficiently large n, a new coloured random graph X with spin law £ and connection probability — , 
such that any edge present in X is also present in X. Let \E\ be the number of edges of X. Using 
Chebyshev's inequality, the binomial formula, and Lemma 12.2} we have that 

n(n — 1) 

'n{n-l)/2\ (C\*(, c\n{n-l)/2-k 

ni 

fe=0 

= e^fl -°- + e ^y {n - 1)/2 < e -n We -i +0 (i))_ 
V n n) 

Now given 9 > choose iV£N such that N > 9 + c(e — 1) and observe that, for sufficiently large n, 

F{\E\ > nN} < F{\E\ > nN} < e~ n0 , 
which implies the statement. ■ 

Lemma 2.4 (Exponential tightness). For every 9 > there exists Kg C 7W(M) such that 

lim sup \ logPJL 1 K e \ < -9. 

Proof. Let 1 < Z G N and choose k(l) G N, large enough, such that £[e l2l i a > k W<} < 2 l 
Then, using exponential Chebyschev's inequality, we have that 

pf / L l {dx) >)] < e- nl E\e l2 ^vhsWu))>H l)} \ < e ~ nl (£[ e l2l {s>m) } ]) n < e -n(J-iog2)_ 

Now we fix 9 > 0, choose M > 9 + log 2, define the set Tm by 



1,—n \ / 



r M := jw : / w(tte) > j, for all Z > Af ) 



'{x>fc(0} 

As {x < k(l)} C 1 is compact, the set Tm is pre-compact in the weak topology, by Prohorov's 
criterion. As 



"{l 1 G* T m } < T -i- 3T exp(-n[M - log 2]), 
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we conclude that 



limsup - logPJL 1 g" K e \ < -6, 
for the closure Kg of Tm as required for the proof. 



2.4 Proof of the upper bound in Theorem 13.41 

We denote by C\ the space of bounded functions on R and by C 2 the space of bounded symmetric 
functions on 1 x R, and define 

I(u,m) = sup { f {f[x\-U f )u[dx]+\ [ g[x,y]w[dx,dy} + \ [ (l-e 9 ^)C[dx,dy]cu[dx]u[dy}\. 

sec 2 
Lemma 2.5. For each closed set F C M(l) x M*(R x R), 



limsup ± log P{(L\L 2 ) EF} < - inf I(u,w). 

n— >oo ' (w,ra)ef 

Proof. First let / € C t and g G C 2 be arbitrary. Define J3 : R x R ->• R by 

/3[x,y] = (l-e^)C[dx,<%]. 

Observe that, by Lemma [272], /3[x,y] = lim n ^. 00 /i n [x,y] for all a, b G R, recalling the definition of /t n 
from (|2,3p . Hence, by (|2.4p . for sufficiently large n, 

e sup xem \p(x,x)\ > [ e (h L A>~ h ") dP = E I e n(L\f-U f -)+n(l J L 2 ,g)+n(±Ll®L\h n )\ 

where L^ = ± E« e v <*(X(«),X(«)) and therefore, 

limsupilogE{e n<Ll ' / - C/ / )+n< ^ 2 ^ +n ^ Ll0Ll ^" ) } < 0. (2.6) 

Given s > let I e (ui, w) = min{I(cj, w),e~ 1 } — e. Suppose that (u, w) G F and observe that I(u, w) > 
I £ (u), w). We now fix / G C\ and <7 G C2 such that 

(/- Upu) + \ (g,xu) + I0,u(g)(jj)> I e (u},w). 

As /,5 are bounded functions, there exist open neighbourhoods B^ and B^ of w,lo such that 

inf {(f-Uf,u>) + \{g,m) + ± (/3,w ® w)} > I e {oj,vj) - e. 

Using Chebyshev's inequality and (j2.6[) we have that 

limsup^ log P{(L\L 2 ) G Si x B^} 
n—^oc 

< hmsup ± logE< e x J /' X2 ' N2 ' > — I £ {co,w) + e (/.7J 

n— >oo *- J 

< — I £ (o;,ro) + e. 

Now we use Lemma 12.31 with 9 = e _1 , to choose N(e) G N and K £ such that 

limsup i log PJ \E\ > nN(e)\ < -e _1 and limsup ^ log PJL 1 K £ \ < -e" 1 (2.8) 
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For this N(e) and K £ we define the set T £ by 

T £ = Uu,m) G M(R) x M*(R x R) : w G iQ c , [|w|| < 2iV(e)}, 

and recall that ||F 2 || = 2\E\/n. The set T e n i 7 is compact and therefore may be covered by finitely 
many sets B^ r x B^ r ,r = 1, . . . , m with (u> r ,zu r ) G F for r = 1, . . . , m. Consequently, 

m 

P^F 2 ) G F} < ]TP{(L\L 2 ) G < x 5^} +P{(L 1 ,L 2 ) ^r e }. 

We may now use (|2.7p and f|2.8j) to obtain, for all sufficiently small e > 0, 

limsupilogPKL 1 ,! 2 ) G F} < max (limsup ± logP{(L 1 , L 2 ) G B^ X S^}) V (-e)" 1 

<(- inf I £ {uj,w) + e]y {-eY 1 . 

V (w,ro)eF / 

Taking e \. we get the desired statement. ■ 

Next, we express the rate function in term of relative entropies, see for example Dembo et al. [1998, 
(2.15)], and consequently show that it is a good rate function. Recall the definition of the function I 
from Theorem 13.41 

Lemma 2.6. 

(i) i{oj,w) = I(cj,w), for any (u,w) G M(M) x M*(R x R) ; 
(ii) I is a good rate function and 
(iii) $jc(w || w) > with equality if and only if w = Cuj (g> to. 

Proof. (i) Suppose that w ^ Ceo (g) ui. Then, there exists ao, &o G R with Cu (g> u(ao, bo) = and 
w(a , b ) > 0. Define g: R x R ->■ R by 

5[x,y] = log [if(l (a0)6o )[x,y] + l {b(hao) [x,y}) + l] , for a, 6 G R and FT > 0. 

For this choice of g and / = we have 

(f[x] - U f )oj{dx) + I \g[x, y}vj[dx, dy] + [ |(1 - e^ ] )C[(ix, dy]u;[dx]w[dy] 

JRxR JExR 

> ilog(i^+ l)tu(ao,6 ) — ^ oo, for FT | oo. 

Now suppose that w <C Cw ® w. We have 

J(w,ot) = sup { / (f[x\ - log / e f[x k[dx]) ui[dx}} 
fed Ur v Jk / J 

C[dx,dy]w[(ix]a;[dy] + \ sup s / g[x,y]w[dx,dy] — / e s ^'^C[(ix,dy]w[dx]a;[<iy] f. 



+ 3 



KWABENA DOKU-AMPONSAH 



By the variational characterization of relative entropy, the first term equals H(ui\\£). By the substi- 
tution h = e 9 ®^ the last term equals 



sup 

hec 2 jRxl 



i ( i r 1 zu[dx,dy] \ i 

l0g { h[x ' v] C\dx.dvMdxWdvV ~ h[xM \ W[dxM 



= — \\vj\\ + H{w || Cw (8 oo), 

where we have used sup x>0 logx — x = — 1 in the last step. This yields that I(co,zu) = I(uj,w). 

(ii) Recall from (j!T2j) that 

I(oj, w) = H(u \\£) + \ H{w (I Cuj w) + I ||Cw 8) w|| - | ||ro||. 

All summands are continuous in to, w and thus / is a rate function. Moreover, for all a < oo, the level 
sets {/ < a} are contained in the bounded set {(uj, w) G M.(M) x 7W*(M x R) : ficfo || w) < «} and 
are therefore compact. Consequently, / is a good rate function. 

(iii) Consider the nonnegative function S,(x) = xlogx — x + 1, for x > 0, £(0) = 1, which has its only 
root in x = 1. Note that 

Sic{w || W ) = { ^° 9dCuJ ® ^ lf 5 := *fe " ° 6XiStS ' (2.9) 

I oo otherwise. 

Hence ficfa II w ) > 0, and, if w = Cuj (g> w, then C(dc^u) = £W = an d so $)c(Cuj <g> w || w) = 0. 
Conversely, if £jc(^||w) = 0, then -cufda^dy] > implies Cw w[x,y] > 0, which then implies 
£ o y[x, y] =0 and further y[x, y] = 1. Hence to = Cuj (g) w, which completes the proof of (iii). ■ 

2.5 Proof of the lower bound in Theorem 13.41 
Lemma 2.7. For every open set O C M(M) x .M*(R x R), 

liminf ± logP^L 1 , L 2 ) G o\ > - inf I(oj,w). 

Proof. Suppose (a;, ru) G O, with ro <C Cw w. Define / w : R — >■ E by 

/ w (a) = { log W' ifw ^]>0, 
I 0, otherwise. 

andfo:lxI^IRby 

I 0, otherwise. 

In addition, we let /3 ro [x,y] = C[dx,dy](l — e 9 ^ x,v ^) and note that /? ro [x,y] = linin^oo /i rojn [x,y], for 
all a, b G R where 

/i ro ,n[x,y] = log (l - p n [dx, dy] + p n [dx,dy]e 9 ™ [x ' y] y n . 

Choose B^,B^ open neighbourhoods o£u,m, such that B* x i? 2 ^ C O and for all (uj,w) G -B^ x fi 2 ^ 

(/w,w) + ^(y ro ,ro) + i(^ w ,w(8)w) -e < (fu,,&) + ^(g^,w) + | (/3 ro ,w 8) w). 
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We now use P, the probability measure obtained by transforming P using the functions f u , g m . Note 
that the spin law in the transformed measure is now u, and the connection probabilities p n [dx,dy] 
satisfy np n [dx,dy] — > w[dx,dy]/(to[dx]uj[dy]) =: C[dx,dy], as n — > oo. Using (|2.4p . we obtain 

p{(L\L 2 ) e o} > E{g (X)t {{ v jL2)eB i xBl} } 

= E{ H e -/«[«^(«))] [] e -3~[«(»j(«)).«foW)] JJ e-«' i -"[*W tt W'«W , '^]l {(i i^. )6Bix ^ } } 

uGV (u,v)eE (u,v)e£ 

>exp(-ra(/L,w)-n±(<? ro ,ro) - n±(/3 ro ,cj ® w) + m - ne) xP^I 1 ,! 2 ) G Bi x £^}, 
where m:=0A inf xe K/3[x,x]. Therefore, by (J2.5J) . we have 
liminf i logPifL 1 ,^) G o) 

> ~(/ w ,w) - § (g w ,w) - \ (~^,u® u)-e + liminf ± logp{(L\ L 2 ) G ^ x B*,}. 

The result follows once we prove that 

liminf I logp{(L\ L 2 ) G ^ x B|) = 0. (2.10) 

We use the upper bound (but now with the law P replaced by P) to prove (|2.10p . Then we obtain 
limsupilogP{(L\L 2 ) G (£* X Blf) < - inf J{Q,w), 

n— >co (iD,ro)eF 

where F = (B* x B^) c and 

I(&, to) := U(u) || w) + 5 H{w || Cw <8> w) + | HCcD (8) w|| — \ \vj\- 

It therefore suffices to show that the infimum is positive. Suppose for contradiction that there exists 
a sequence (ui n , w n ) G F with I(u n , vb n ) I 0. Then, because I is a good rate function and its level sets 
are compact, and by lower semicontinuity of the mapping (a), tu) t-t I(tD, w), we can construct a limit 
point (Co,vj) G F with I(u),xb) = . By Lemma 12.61 this implies H(Q \\u) = and Sjci^ \\&) = 0, 
hence ui = cu, and zu = CCj ®6j = zo contradicting (a), zu) E F. ■ 



3. Results and Discussion. 



3.1 Thermodynamic Limits 



Theorem 3.1. Suppose that X is a spinned random graph with bond probabilities p n : X x X — > 
[0, 1] satisfying np n [s(x), s(y)] — »■ C[s(x), s(y)], /or a// x,y £ {—1, 1} a^rf some continuous symmetric 
function C:Kxl4 [0, oo) i/iai satisfies (|2.ip above. Then, the asymptotic pressure per particle is 



<P(P, B) = \[{e? - l)C[y/pB(l), y/pB(\)\ + (e^ - 1)C[-^B(-1), y/0B(l)]} 



(3B(1)-B(-1)) 

"r 4 



Corollary 3.2. Suppose that X is a spinned random graph with bond probabilities p n : X x A' — > [0, 1] 
satisfying np n [s(x),s(y)] — > C[s(x),s(y)], for all x,y G {—1, 1} and some twice differential symmetric 
function C:lxl-> [0, oo) that satisfies (I2.1D above. Then, for all (3 > 0, eac/i o/ t/ie following 
statements holds a.s.: 
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(i) Asymptotic total magnetization per site is given by 



M(/3, B) = ~- — ' "V : *n , —^— 

2 ^(^ ) i-^-i,-i)+(^-i)(% i -^r i ) 

(ii) Asymptotic internal energy is given by 

U(P,B) = -\ [C-^eP - C-^e-P + (e-/> - 1)^ + (e' 3 - l)^pi 
(iii) Asymptotic specific heat is given by 

Remark 1 : Asymptotic susceptibility of the system is given by 



&b(-i) 2 as(i) 2 

Theorem 3.3. Le£ /3 C 5e i/ie solution of equation C-\ t \ = \C\^ + ^C-x,-i- Then, (3 = j3 c , separates 
the non magnetized phase fdiamagnetic state,) from the magnetized one ^paramagnetic state). Further, 
for -B(l) + B{— 1) = ; the continuous function B — > <p(f3,B) exhibits a discontinuous derivative. i.e. 
the susceptibility becomes infinite at (3 = (3 C in B{\) + B{— 1) = 0. 

Corollary 3.4. Let A be the asymptotic average number of connectivity of erdos-renyi graphs. Then, 
there exists f3 c (\) such the ferromagnetic Ising model on the erdos-renyi graphs exhibit 

(i) diamagnetic properties if < (3 < /3 C (A) 
(ii) paramagnetic properties if (3 C (\) < (3 < oo 
(iii) zero-magnetization properties if j3 = or (3 = f3 c (X). 

3.2 Proof of Thermodynamics Limits. We kick start the proof of our Thermodynamics limit 
results by stating an important Lemma (Varadhan's Lemma, see Dembo et al. [1998 „ Theorem 4.3.1] 
) , which is a key step in establishing our first result Theorem 13.11 without proof. 

Lemma 3.5 (Varadhan). Suppose the functional M n from the space of finite spinned graphs to M 
satisfies the LDP withh good rate function I : M. — >• [0, oo] and let *$> : Ai — > R be any continuous 
function. Assume further the following moment condition for some A > 1, 



lim sup - log E [e 



n\V{M n [X])] 



< OO. 



Then, 

lim IlogEre nA * {M "™l = sup fV(m) - I(m)\. 

Next we provide annealed asymptotics of the random partition fuction for the ferromagnetic Ising 
model on spinned random graphs, as the graph size goes to infinity. 

Lemma 3.6. Suppose that X is a spinned random graph with spin law £ and link probabilities 
p n : X x X — > [0,1] satisfying np n [s(a),s(b)] — > C[s(a),s(b)], for all a,b S {—1,1} and some con- 
tinuous symmetric function C: R x R — > [0,oo) with sup xy ^C[x, y] < oo and satisfies (|2.ip above. 
Then, 

lim I log E[Z x (f3,B)] = log 2 + \{{e? - 1)C[^B(1), y/pB(l)} + (e~ p - 1)C[-^B(-1), y/0B(l)]} 

(3B(1)-B(-1)) 



ISING MODELS ON STATIC INHOMOGENEOUS RANDOM GRAPHS 11 

Proof. Recall the random partition function of the Ising model on the spinned graph from section [5] 
and write it as integral of some function with respect to our empirical measures: 

E[Zx(/3, B)] := 2"E [exp {§ J ^^ L 2 [ds(x),d S (y)] +nj ^L l [ds(x)]}\ (3.1) 

Now using the (Varadhan) Lemma 13.51 and Theorem 13.41 we obtain, 

Jim I log E[Z x (/3, B)] = log 2 + sup [\ J ^^ w [ds{x), ds{y)\ + J s -§u[ds(x)} - I(u, w) 

: u G X({s(-l),s(l)}),w G M*({s(-l),s(l)} x { a (-l), s (l)})} 
= sup|f (w(A) - zu(A c )) + B(l)x - B(-l)(l - x) -xlog(x) - (1 - x) log(l - x) 

- § (#(ro || wj.) + Ci,ix + C_i _i(l - x) + 2C_i,ix(l - x) - ||w||) } , 

(3.2) 
where A is the diagonal in {s(— l),s(l)} x {s(— l),s(l)}, and the supremum is over all x G [0, 1] and 
w G M*({s(— l),s(l)} x {s(— l),s(l)}), and the measure ui x G .M*({s(-l),s(l)} x {s(— l),s(l)}) is 
defined by 

u s [s(i),8(j)] = Cyi^+^a-x)^-^ 2 for i,j G {-1,1}. 

We take the partial derivatives of (13.21) with respect to w, and set our results to zero to get, 

ro [ S (l) lS (l)]=eV[s(l),s(l)], 

zu[s(-l),s(-l)} = e" Wl [ S (-l), S (-l)] 

tJ7[a(l),a(-l)]=e-^[ S (l),s(-l)] 

w[s(-l),s(l)}=e-Poj x [s(-l),s(l)]. 

Hm ±logE[Z x (/?,5)] 

= sup { f A; x [s(l), a(l)] + | eV[s(-l), a(-l) - /^^(-l), s(l)] + S(l)x - S(-l)(l - x) 
xe[o,i] L 

- xlog(x) - (1 - x) log(l - x) - | A> x [fl(l), a(l)] - | e^ x [s(-l), s(-l)] + /Se^^-l), a(l)] - ±C\ix 

- |C_i,_i(l - x) - C_ 1)1 x(l - x) + ieV[s(l)> *(1)] + l^K-1), s(-l)] + e-%[s(-l), s(l)]} 

= sup { -xlog(x)- (1 - x) log(l - x) - ±Ci,ix 2 -±C_i_i(l-x) 2 + £(l)x-5(-l)(l-x) 
ie[o,i] *• 

- C_i,ix(l - x) + \e p C lt ix 2 + |e^C_i,_i(l - x) 2 + e^C-i,^! - x)} 



= sup I Ai(x) + A 2 (x) > 
ase[o,i] L J 

where 

Ai(x) = — xlog(x) — (1 — x) log(l — x) 
A 2 (x) = aix 2 + o 2 (l - a;) 2 + a 3 x(l - x) + B(l)x - S(-l)(l - x) 
ai = ± Ci,i(e" - 1), a 2 = § C_i,_i(e^ - 1) and a 3 = CLi,i(e-* - 1). 



(3.3) 
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Elementary calculus shows that the grobal maxima of the functions Ai and A2 are attained at the 
value x = 5, when 02 — a\ = B(l) + B(— 1). Therefore we have 

lim ± log E[Z X (/?,£)] = log 2+ i(2oi +a 3 ) + i(3S(l) - fl(-l)) 



± |V - 1)C[V#B(1), V^B(l)] + (e-f 3 - 1)C[-^B(-1), V^B(l)]} 



Lemma 3.7. 



$ n Cg,i?)-log2 MogE[Z x (/3,i?)/2"] 
hm 2 = urn j with prob. 1 



Proof. 

^ $ n ((3,B)-log2 ^logE[Z x (/3,B)/2 n ] ^ ; B mB y log2+£ A 

In* : — \ -; — y -+- — it -+* T \ — :r 

I n i n i I — An 1 n A n 4 — n 2 

for some A = A(/3, -B, e). 

Now 



■X 



f ,, «,(/». b) - io g 2 _ 1 log g& eg. gvg , <T 4<oo 



n 

ra=l n=l 



Using the Markov's inequality and ()3,7p we have that, 



OC 



$ n (^,5)-log2 ±logE[Z x (fti5)/2»]. -, .^ ,$„(£, B)- log 2 ±logE[Z x ((3,B)/2 n ] 



Z^ \\ n 4 n 4 I - J - s 2-> I n 4 

n=l n=l 



E 



00 



<iV4<co 

n=l 

By the Borel-Cantelli Lemma, this estimate result in 

$ n ((3,B)-\og2 ilogE[Z x (/?,ff)/2"] 
lim 7 = lim j almost surely 

n— >oo n n— >oo tt 4 

We use Lemmas (|3.2p and (|3.7p . to obtain 

lim *„(A*) = I- n 4^,B)-k>g2 + log2 = Hm ^ ^ * M B) ,- log 2 + 2 

4 ilogE[Zx(^B)/2-] 
hm n hm ( 2 ) + log 2 



n— »oo n^-oo 



= hm ilogE[Zx(/3,i?)]) 

n— »oo 

= <t>{P,B) 
with probability 1, which concludes the proof of our main Thermodynamic Limiting results. 
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Corollary I3T21 

(i) We take the partial derivative of (j>(/3,B), using the chain rule of differentiation, with respect to 
B(l) and B(-l) to get 

alfe = \ [tfC^ - e^Cu) + {e?-l)^ + (e^ - l)*fr] ^ - \ (3.5) 

and 

a^iy = \ [&C ltl - e-'C- ltl ) + (** - 1)-^ + (e-P - 1)«%±] ^ + |. (3.6) 

Now the energetic preference condition gives 

afry = a#T) = -2^1,1 - CU,-i) + (^ - 1)(^ - ^)] _1 . (3.7) 

We write this expression agrn in (|3.5|) and (|3.6|) . and adding them together we have the M(/3, £>). 

(ii) We take the partial derivative of (f)((3, B), using the chain rule of differentiation, with respect to f3 
to get 

I = \ [&Ci,x ~ ^-1,1) + {ef> - l)*fr + (e^ - 1)^] J^ + i(3^ - ^). (3.8) 

Now writing (|3.7p in (|3.8|) and negating our results we obtain U(/3,B). 

(iii) Differentiating U with respect to /3 and write the expression in the relation 



w,*> = -f% 



we have the asymptotic specific heat. 



Proof of Theorem 13. 3t Let /3 C be the solution of equation C-\^_ = gCi.i + op—\,—\- Then, it is not 
uneasy to see that (M(/? c , B) = 0. Further, B(l) + S(— 1) = 0, implies C\ t \ = C-\-\ by the energetic 
preference condition, and /3 — > M(/3,B), the derivative of 4> with respect to /3 is discontinuous for all 
B. Since the susceptibility is the derivative of M(/3, B) with respect to /3, we have 

dM 
hm — — - = oo. 

p^p c d/3 



Proof of Corollary !3.4t Let A be the asymptotic average number of connectivity of erdos-renyi 
graphs. Then, we have np n [s(x),s(y)\ — > A, for all x, y{— 1, 1}. Now take C\^\ = C—\ t \ = C-i-i, and 
note that C-\ t \ = \C\^\ + ^CL^-i, for all (3 > 0. Using Theorem 13.31 we can find /3 C (A) with separates 
the paramagnetic state from the diamagnetic state which ends the proof of the corollary. 
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